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ABSTRACT 

In this dissertation, we consider the weighted goal 
programming and interval goal programming problems in the 
network context. Special structure embedded in their problems 
is used to develop computationally efficient algorithms. In 
weighted goal network flow [WGNF] problem, the decision 
maker specifies aspiration levels for each of the objectives 
and weighting factors for each of the deviations from the 
aspiration levels. On the other hand, in interval goal network 
flow [IGNF] problem, the decision maker provides a range of 
aspiration levels for each objective and weighting factors 
for each of the deviations from the specified range of aspi- 
ration levels. Both the algorithms developed trace a path 
in the feasible region and obtain the solutions which minimize 
the sum of weighted deviations. The WGNF and the IGNF algori- 
thms are based on the parametric approach and fully exploit 
the special structure of the minimum cost flow problem in 
order to perform all the computations on the network itself. 
Computer programs were written for both the algorithms, and 
tested on randomly generated network problems. Results of 
their investigations are presented. 



CHAPTER I 


INTRODUCTION 


I ntroduc tions 

Operations Research is a scientific approach to solve 
complex problems arising in the management of large systems of 
men, machines, materials and money. Operations Research is a 
decision science which helps management to make better decisions. 
Network model .Is a branch of study in Operations Research, Due 
to the wide applicability of network models in real world, it 
is considered as an important branch of study in Operations 
Research. In the present day, we find that complex intriguing 
problems arising in production-distribution systems, military 
logistics systems, urban traffic systems, railway systems, 
communication systems, pipeline network systems, facilities 
location systems, file merge systems, electrical networks etc. 
can be tackled by constructing an appropriate network model. 
Furthermore, network geometry (or relationships) can be easily 
displayed in two-dimensional drawings, greatly simplifying the 
communication problem between the analyst and the client for 
whom the model is designed. In the recent years, network flow 
problems have received special attention due to significant 
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advances in implementation technology and solution techniques, 
thereby increasing the applicability of the network models 
substantially. 

In this dissertation, we consider the bicriteria 
minimum cost flow problems and propose two algorithms i,e, 
weighted goal network flow (WGNF) and interval goal network 
flow (IGNF) algorithms for obtaining optimum solutions. These 
models are useful when two conflicting objectives are to be 
simultaneously considered, i.e, one is interested in minimiz- 
ing the total cost as well as total time. 

One application of the bicriteria minimum cost flow 
problem arises in a distribution system. Suppose V units of 
perishable items (i.e, fruits) are to be sent from source to 
sink along the road network. The two objectives to be considered 
simultaneously are total cost and total time. We would like to 
transport the items in such a way that they are not perished 
before they reach the sink and at the same time the total cost 
of transportation should be minimum. 

Another application of the bicriteria minimum cost flow 
problem arises in communication networks. Suppose V units of 
messages/sec. are required to be transmitted from source to 
sink along a communication network. The decision maker would 
like to consider two objectives simultaneously. One of the 
objectives may be to minimize total cost of sending messages, 
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and the other being the minimizat l r n of total amplification 
cost. 

Another important application of bicriteria network 
flow problem is in the area of flow through a pipeline network. 

The decision maker would like to minimize the total operating 
cost for sending the flow as well as minimize the total time 
to transmit flow from source to sink. 

In this work, we develop specializations of goal 
programming i.e. weighted goal and interval goal network flow 
algorithms based on parametric approach. The special structure 
of minimum cost flow problem is exploited to perform all the 
computations on the network itself. Efficient tree data structures 
are used to further enhance the efficiency of these computations. 
Computational investigations with the above approaches are found 
to be very encouraging and presented in sufficient detail. 

1.2 Outl ine, of the^ .Thesis; 

In this section, we give a brief outline of this 

thesis. 

In Chapter II, we present some of the preliminary 
knowledge required to understand the work done in the subse- 
quent chapters. We first review the literature related to the 
bicriteria minimum cost flow problem. The graph theory nota- 
tions adopted in this work are also described. Finally, an 
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overview of the bicriteria nebwcrK flow problem and the goal 
programming are briefly outlined. 

Chapter III embodies the main work done on WGNF problem 
and IGNF problem. They can be formulated as shown belows 


Weighted Goal Network Flow Problem^ 


Min. Z = w 2. a l + W 2 a 2 + r l^l * r 2^2 

s.t. 


. -V, if i=] 

\ 9 


Z x.. - Z x. = ; V, if i=n, 

(j,i)el(i) 31 (i,d)eOCi) 10 / ¥ xeN 

V 0, Otherwi 


0 < x ij $ Vi* V (i»o) e A 


Z c . x. p, - a n = C 

a, aw- 13 13 1 1 


2 d. • x* » *■** (3^ *•* 
(i,j)eA 13 13 2 2 


Interval Goal Network Flow Problem ° 


Min. Z = w^ -i- w^ a ^ + r 1 P 1 + r 2 P 2 


s.t. 


f-V, if i = 1 

Z x. . =s j V, if i = n, V iej 

(j,x)el(i) (i,o)eO(i) ^ q, otherwise 


X . . 


0 < x . .. < b _ , V (i,o) e A 

ij ij ’ 


Z c . x* *1" p-i — oc-i — C-. 

(i,a)eA 13 13 1 1 1 


? . d ,-, x <, + Po - “o - D i 


(i,d)eA ^ ^ H 2 2 
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Z c x. 
(i,Ct)eA 10 10 


a- 


= C, 


I 

(i,o)eA 






where, 

x. . is the amount of flow on arc (i,j). 

XJ 

c. is the cost of unit flow on arc (i,j), 

X J 

d. . is the budget required for unit flow on arc (1,3) 

b^ is the capacity of arc (1,3). 

V is the net flow of commodity from source to 

sink n. 

Ch is the aspiration level for cost ob3ective function. 

D i is the aspiration level for budget ob3ective 

function. 

is the positive deviation, 
is the negative deviation. 

We suggest two exact algorithms to solve bicriteria network 
flow problem. They are weighted goal network flow (WGNF) and 
interval goal network flow (IGNF) algorithms. Both these 
algorithms are based on the parametric algorithm for constrained 
minimum cost flow problem [l] . 


The WGNF algorithm searches for an optimum solution 
which satisfies the aspiration levels specified for each objec- 
tive by the decision maker. Whereas IGNF algorithm searches for 
an optima that lies in the specified range of aspiration levels 
for each of the objectives. 
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The computational investigations of both the algorithms 
are given at the end of Chapter III. Computer programs are 
■written for both the algorithms and the computational perfor- 
mances are ascertained by solving different sized problems. 
Efficient data structures are implemented to represent the basis 
tree which requires comparatively less storage. Data struc- 
tures based on augumented Threaded index method are used and 
results of the computational investigations are presented* 



CHAPTER I" 


PRELIMINARIES 


2.1 Introd uct ion s 

In this chapter, we review some of the relevent concepts 
of the Licriteria minimum cost flow problem and goal programming. 
The basic reason behind this is to prepare a sufficient back- 
ground which will help in understanding the algorithms developed 
m the subsequent chapters. 

This chapter is divided into five sections, The graph 
theory notations adopted in this work are presented in Sec. 2. 

The literature related to the bicriteria minimum cost flow- 
problem is surveyed m Sec, 3. A brief review of various 
methods of solving bicriteria minimum cost flow problem and a 
parametric algorithm for solving constrained minimum cost flow 
problem [l] are presented in Sec. 4. In Sec, 5 a brief des- 
cription about goal programming and its importance to manage- 
ment in decision making to real world situations, and various 
goal programming techniques are discussed. 

2.2 Graph The ory Not atipns '* 

Some notations and well-known concepts of graph theory 
that are being used throughout the thesis are mentioned below. 
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A dlrected^jragh G = (1 m s A), consists of a finite 
set N of elements* called no des., and a set A of ordered pairs 
of nodes called arcs,. A directed, network is a directed graph 
in which numerical values are attached to the nodes and arcs 
of the graph* Let n = |n| and m = 1 A j ^ The two specified nodes 
1 and n are called the sourpe, and the sink respectively. 

An arc (i,o) has two end points, i and j, and it is 
said to be incident, from, node i and incident^to. node j . Let 
l(i) and 0(i) denote, respectively, the sets of arcs incident 
to and incident from node l. The degree of a node i is the 
number of arcs incident to or incident from that node. 

A p ath in G = (N, A) is a sequence i^, i 2 , ...» i r 
of distinct nodes of N such that either (± k , i k+ ^) e A or 
(ifc+1* i]P e A for each k = 1 ,..,, r-1 . A directed^path is 
defined similarily, except that (i k , i k+ -^) £ A i° r eac ^ 
k = 1,.*., r-l. A cycle is a path together with an arc 
(i r , i x ) or (i^, i r ). A directed, cycle, is a directed path 
together with the arc (i r , i^) . 

A graph G = (N' # A' ) is a subgraph of G = (N, A) if 
N'c N and A' qA, A graph G » (N' , A' ) is a spanning subgraph 
of G = (N, A) if N' = N and A’c A # 

Two nodes i and j are said to be connected if there 
is atleast one path between them. A graph is said to be 


t 
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connected, if all pairs of nodes a >e connected, otherwise it 
is called disconnected. A set Qc A such that the graph 
G s= (N, A-G) disconnected and no subset of Q has this property, 
is called a cocycle of G, A cocycle is a cutset if it dis- 
connects source and sink, 

A graph is acyclic if it does not contain any cycle . 

A tree is a connected acyclic graph, A subtree of a tree T 
is a subgraph of T as well as a tree, A tree T is said to be 
a scann ing tre e of G if T is a spanning subgraph of G. Arcs 
belonging to a spanning tree T are called tree-ar^cs, and arcs 
not belonging to T are called nontree-arc s_ c A spanning tree 
of G s (N, A) has exactly (n-l) tiee-arcs. 

2 . 3 Literature Review « 

One of the methods of solving the bicriteria minimum 
cost flow problem is to consider one of the objectives as a cons- 
traint and solve the resultant constrained minimum cost flow 
[CMCF] problem. We will first present a brief review of the 
constrained flow problems. 

Hultz and Klingman [9j have suggested a partitioning 
method in conjunction with the simplex method for solving 
constrained generalized network flow problems. Takashi 
Kobayashi [l2] proposed a primal-dual method for solving CMCF 
problem. He has associated two dimensional distance for each 
arc. The first element is related to the cost and the second 
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one to the coefficient of the ado 5 -tonal constraint. This 
method is suitable m cases when degeneracies often occur, 

A parametric algorithm for solving CMCF problem is 
developed by Ahuja, Batra and Gupta [l] , The algorithm uses 
the concepts from parametric linear programming and fully 
exploits the topological scructure embedded in the problem. 

The algorithm can also be used to generate all the efficient 
points of the bicriteria minimum cost flow problem. 

Chen and Saigal [4] have suggested a primal algorithm 
for solving a capacitated network flow problems with additional 
linear constraints. 

Klmgman and Mote [10] have reviewed the fundamental 
theoretical results for the general multi-criteria linear 
programming and the relevent exploitable characteristics of 
the network basis. Two algorithms are then developed for solv- 
ing the multicriteria network flow problems efficiently. One 
approach determines the set of all non-dommated solutions to 
the problem. The other approach is a network variant of the 
surrogate criterion linear programming approach. Ignizio [7] 
has suggested a straight forward weighted integer goal programm- 
ing for generalized networks models for integer programming 
problems. His method is simplex and robust. 
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Practically no study wap lacle in + t^ field of appli- 
cation of goal programming techniques to network flow problems. 
This has motivated us to propose two algorithms for WJNF and 
IGNF problems. 

2.4 Bicriteria Line ar P rogramming; 

Bicriteria linear programming deals with the optimiza- 
tion ol two objective functions (multicnteria deals with two 
or more than two ^objective functions) simultaneously. A deci- 
sion situation is generally characterised by multiple objectives. 
Some of these objectives may be complementary, while others may 
be conflicting in nature. 


The bicriteria minimum cost flow problem differs from 
the classical minimum cost flow problem only in the expression 
of their objective functions. The bicriteria minimum cost flow 
problem is shown as belowi 


Minimize 



subject to 

£ x . 
( j, i)el(i) 31 


£ x. 
(i, j)eO(i) 



if i=l 

if i=n, V ieN 
otherwise 


0 ^- x i0 s-V 


v (i,j) e A 
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Generally, it is observe] that these two objectives 
are conflicting in nature and hence there is no optimal solution 
in the normal sense for thd above problem* The decision maker 
has to choose solutions* possibly not the best for both the 
criteria. A special set of solutions, the non- dominat ed or 
ef fief e_nt_ solutions can be defined to overcome this problem* 

An e f f i ciert _so ly cion is one in which one objective 
cannot be reduced without a simultaneous increment in the other 
objective* That is, X* is an efficient solution to the bicri- 
terion minimum cost flow problem if there does not exist any 
X° e S, the set of all possible solutions such that 

Z-^X 0 ) £ Z^X* ) and 
z 2 (x°) < Z 2 (X* ) 

with atleast one strict inequality. 

A number of techniques as mentioned below are available 
for generating efficient or non-inferior solutions for the 
bicriteria problem formulated above. 

(i) The weighting method 

(ii-) The constraint method 

(iii) The non- inferior set estimation method, and 

(iv) The multi-objective simplex algorithm. 

These techniques are discussed in detail by Ambrose Goicoechea [3] . 
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There is one disadvantage with non- nominated solutions 
technique. It becomes difficult for the decision maker to 
make his final choice from 3 set of non- dominated solutions. 
Another technique known as goal programming, allows the deci- 
sion maker to specify a target for each of the objective 
functions. It obtains a p >jref errpd,_splutipn which is defined 
as the one that minimizes vhe sum of the deviations from the 
prescribed set of target values. A brief description of the 
goal programming is presented m Sec. 2.5. 

An_ pveryi_ew_ of_ Co ns trained^ Minimum, Cost Flow Problem* 

We now briefly outline tie parametric algorithm for 
the constrained minimum cost flow [CMCF] problem which is the 
basis for the algorithms developed in this thesis. 


The mathematical formulation of the CMCF problem is 


Min. 


s.t. 


Z = £ c 

(1,3 )eA 

£ x. . - 
(3, l)sl(i) 


0 x l 3 b ic’ 

£ d. x. 
(i,j)eA 10 10 


ij 

1 x . 
( 1 , j)eO(i) ^ 

¥ (i, j) e A 

< D 


( 2 . 1 ) 


/-V,if i=l 

J V,if i=n, V ieN 

/ 0, otherwise 

1 ( 2 . 2 ) 

(2.3) 

(2.4) 


This algorithm is developed by Ahuja, Batra and Gupta [l]. The 
algorithm utilizes the concepts of parametric linear programming 
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and performs all the computation? over the network. It first 
obtains the optimum basis corre so ending to where 

is the minimum value of C for which a feasible flow exists 
and then moves from one optimum basis to the next as D is 
decreased. The algorithm thus yields the minimum cost as a 
function of the budget, which is a piecewise linear convex 
function* 

Optimality Conditions 5 

It is Veil known that the optimum basis of MCF problem 

is a spanning tree [ 5 ] * The addition of the budget constraint 

introduces one more arc in the optimum basis. It follows from 

bounded variable linear programming that the necessary and 

sufficient conditions for a feasible basis structure to be an 

optimum basis structure are that the.re exist dual variables 

tc.'s and ^ 0, satisfying the following conditions* 

J 


(i) 

n - 

0 

*i 

{( 

0 

H* 

CJ 

+ 

v U,d) 

s B 

(2.5) 

(ii) 

n - 

0 

*i 

~ C i 0 

+ ^ d i0- 

H* 

CJ 

£ L 

(2.6) 

(iii) 

% - 

0 

n l 

^ C ij 

+ it d . . , 

-LJ 

V U, 0 ) 

£ U 

(2.7) 

(iv) 

A (D 


Z ( 

i x. ) 1 

= 0 


(2.8) 

( 

i,j)eA 

1 j io 




where B, L, and U represent the sets of arcs corresponding to 
the basic variables, the nonbasic variables at their lover 
bounds, and the nonbasic variables at their upper bounds, 
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respectively. The set B is called a basis and the triple 
(B, L» U) is called a basis structure. 

Let (T u t(p,q)}, L, U) be an optimum basis structure 

at D a D 1 , Further let tc c and n d be the numbers satisfying 

J J 

\ C = 0 and - % ± c = c ±y ¥ (i*j) e T (2,9) 

= '0 and 7t 0 d ~ = d ±0 , ¥ (i, 3 ) e T (2,10) 

Define the numbers c\.. and cL.. for all arcs (i,j) e L U U as 
follows 2 




-*d° 

+ c - . 

if 

(i,i) 

- 7C.° 
1 

* c ia * 

if 

(i»o) 

„ d 

+ d id- 

if 

(i»l) 

- *. d 
1 

■ d ij' 

if 

(i.o) 


e L 
e U 

e L 
e U 


( 2 . 11 ) 


( 2 , 12 ) 


Clearly, c. = d = 0, for all (i,j) e T. If be is a real 

X J X J 

number and n be the numbers satisfying, 

u 

^ = 0 and itj - 7t i = c^ -wd.^, ¥ (i,j) e T 

then, 

rt. = 7t c + d ¥ j e N (2.1?) 

J J J 


substituting (2.13) in the conditions (2.5) - (2.8) and then 
using (2,ll) and (2.12), we obtain the following equivalent 
conditions 2 
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(i) 

c" + Ak d =0 

(2.14) 

pq. pq 

(ii) 

c i0 +J( d\ > 0. V (i,j) £ L U U 

(2.15) 

(iii) 

Xl(D' - E d. x ) = 0 

(2.16) 

(i,0)eA ^ ^ 


The conditions (2.14) - (2,16) are subsequently referred to as 
the optimality-^onditj-ons^ The c. . . , 3T . and )J are referred to 
as the cost price, budget price and price ratio, respectively, 
associated with the basis structure (T U i(p,q)}, L, U) . 

Depending upon the values of cl. and d. , arcs in L U U 

-**U X U 

can be classified as follows s 

(i) active arcs? e L LI U: c. . > 0 and d. < 0}, 

-4 0 X v 

(ii) critical arcs: i. (1,0) £ L U Us c 0 and d. s 0}. 

XJ "■“* XJ 

(lii) passive arcs* l(i,o) £ L U U: d 0} , 

X J 

The piecewise linear convex curve between cost and 
budget is obtained as follows. 

Firstly an arc (p, q) belonging to the critical set is 
selected and added to basis. This results in the formation 
of exactly one cycle. Next the flow is augumented in this 
cycle. One of the arcs m the cycles reaches one of its 
bounds and leaves the basis respectively. The dual variables 
as well as (B, L, U) are updated and the step is repeated 
until the critical set is empty. We now describe how the 
algorithm moves from one basis to the next basis satisfying 
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the optimality conditions vlult decreasing the total budget D. 


-Addition of an arc (i,qy £ L U U to the basic tree 
creates exactly one cycle W consisting of the basic arcs. 

We define the orientation of the cycle V. along (i,o) if 

J- J 

(i,j) e L and opposite to (i,j) if (i,d) £ U. Let W. and 
W be the sets of arcs in the cycle along and opposite 
to its orientation, respectively. Then using (2.9) and (2.10) 
it can be easily shown that 


°id = 


2 c. - 2 c . 

(i,o)eW ia (i ,o)eW 


(2.17) 


10 


10 


10 


2 d . - 2 d . 

(i,o)etf ia 0 (i,j)eW aJ 10 


(2.18) 


Thus, c. (or d. ) denotes the increase in the cost (or budget) 

1 J 1 J 

if unit amount of additional flow is circulated in the cycle 
W. . along its orientation. The above given classification of 

X J 

non-tree arcs can be given the physical interpretation. Active 
arcs are those arcs which can lead to increase in cost if 
budget is decreased. Critical arcs are those which do not 
increase cost if budget is reduced. Passive arcs do not lead 
to decrease m budget even if more cost occurs. 


Char act eri stic Interval} 

Let at D = D, the optimum lasts structure of the CMCF 
problem is (T U [(p,q)}, L, U). Let x. . be the flow on arc 
(i»o) £ A and Z be the cost of this flow. Further, let 
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» d • . and be the cost prices; budget prices and price 
ratio associated with the given oasis structure* 

We now determine lie interval (D, EQ for the values 
of D for which the given basis structure continues to remain 
optimum* This interval is known as the char act eristic inter- 
val associated with (T 6 {(p,q)}» L, U)* 

Since the numbers, c d. and are uniquely deter- 

X J X J 

mined for a given basis structure, the optimality conditions 
(2.14) and (2*15) are not affected by decrease in the value of D. 
However, since AL > 0, the flow must be changed in order to 
satisfy (2.16), The only way to change the flow, without chang- 
ing the basis structure and satisfying the flow conservation 
constraints (2*2), is by circulating some additional flow in 

the cycle W along its orientation. It was noted that d <0 
P4 pq 

is the rate at which the budget is reduced and c > 0 is the 

pq 

rate at which the additional cost is incurred when unit amount 
of flow is circulated. Since the changed flow must also satisfy 
the bound restrictions of the arcs (2.3), we calculate the 
maximum flow f that can be circulated without violating the 
bound restrictions of the arcs in W . If we define, 

Jr M. 





if U , 3 ) e 

if (i,j) c W pq 
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then, 


f = min 




pa 


Thus, 


£) = D + f d 


pa 


For all values of D in (D, IT), the optimum flow is given by. 


X i0 J ' a 1 ’ 

x ia “ a ’ f ’ 

x ia 


if (I,j) £ Vj 

if (i,a) £ W 


pa 


otherwise 


■pa 


and. cost by 

Z = 

where, a’ = 


- + F °pq’ 

(D - D)/(tT - D) 


Let x. . denote the flow in arc (i,j) e A at D = £>, If it is 
required to find the optimum flow for D < E>, then a dual simplex 
iteration is performed to obtain a new basis structure at D = D. 


The^ Dual^ Simplex^ ^Iteration* 

At D = If, flow in an arc (u, v) e W , for which 
f* uv as f, equals its lower or upper bound. If D is decreased 
further, the bound is violated. Thus, to obtain a new basis 
structure at D = C, which may permit decrease in the value of 
D, the arc (u, v) is dropped from the basis and a non— basic arc 
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is selected to enter the basis. T] ^ arc (u. v) becomes a 
non-basic arc at its respecxive bound. In the new basis 
structure, the dashed values and sets represent the corres- 
ponding values and sets of the previous basis structure. By- 
selecting an active arc and performing dual simplex iteration 
each time until the active set is empty, the cost is increased 
and the budget is reduced. 

The curve starts at C _ and as the value of D is 

mm 

decreased, the value of Z keeps increasing until a value of 
D is reached when Z stops increasing. The slope of the curve 
at any point is the value of - LL in the optimum solution which 
corresponds to that point. Since the value of -it keeps 
increasing and finally becomes zero and no critical arcs are 
formed at any iteration (The proof is given in [l}), the curve 
is a piecewise linear convex function. 

2*5 An_ Overview of Goal Programming » 

As mentioned in the previous section, Goal programming 
allows the decision maker to specify a target for each objective 
function which provides him the preferred solution by minimiz- 
ing the sum of the deviations from the prescribed set of target 
values . 

Generalized goal programming has a number of special 
terms and concepts that are being used in this thesis. They 
are mentioned as below! 
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(i) Objectives 

An objective is a relatively general statement (in 
narrative or quantitative cerms) that reflects the desires 
of the decision maker. For example, one may wish to ’’maximize 
profit" or "minimize total wastage" or "wipe out poverty ", 

( ii ) Aspiration^ Level s 

An aspiration level is a specific value associated 
with the desired or acceptable level of achievement of an 
objective. Thus, an aspiration level is used to measure, the 
achievement of an objective and generally serves to ''anchor" 
the objective to reality. 

(ill) Goals 

An objective m conjunction with an aspiration level 
is termed a goal. For example, we may wish to "achieve at 
least X units of profits" or "reduce the rate of inflation 
by Y percent." 

(iv) Goal Deviations 

The difference between what one accomplishes and 
what one aspires to is the deviation from his goal. A devia- 
tion can represent over, as well as under\ achievement, of a 
goal. 
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(v) Goal Fo rmulations 

We will now examine how to mathematically transform 
an objective into a goal within our goal programming frame 
work* 

Consider the objective function expressed, in general 
terms as f^(X). The procedure to be presented is applicable 
whether f^(X) is linear or nonlinear, but only linear objectives 
are considered in this thesis. 

We then let, 

f^(X) = mathematical representation of objective i 

as a function of the decision variables 
X = (x-p x 2 , . . x^) 

b^ = value of the aspiration level associated 
with objective i c 

Three possible forms of goals may then result! 

(i) f (X) <_ b^s that is, we wish to have a value of 
f ^ ( X) that is equal to or less than b^. 

(ii) f^(X) X that is, we wish to have a value of 

f^(X) that is equal to or greater than b^, 

(iii) f. (X) = b % that is f (X) must exactly equal b.. 

X X Jl X 

Regardless of the form, we shall transform any of 
these relations into the goal programming format by adding 
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a negative deviation variable ( (3j 0) and subtracting a 
a positive deviation variable (a n - >, 0). This statement is 
summarized as belows 


Goal Type 


Goal Programming Form 


Deviation Variables 
to be Minimized 


f,(X) <. 

= b ± 


f ± (x) - p x 

f,(x) + p. 


f ± (x) + P 1 


a . 


= b. 


a 


- a* - b 


- a. 




1 

h 


“l + h 


We will briefly discuss the three extensions of goal 
programming, They are* 

(1) Weighted goal programming. 

(2) Interval goal programming. 

(3) Fuzzy goal programming. 

W eight ed Goal Programming 

In this model the decision maker assigns an aspiration 
level for each of the objectives and also the weighting factors 
for each of the deviations. It obtains an optimum solution by 
minimizing the sum of weighted deviations. 

There is an alternative way in which the weighted 


model may be formed. Rather than multiplying each deviation 
variable by a constant weight, we may instead, raise each 
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deviation variable m the acnieve lent function to some power. 
This results in a polynomial form for the achievement function. 

Given a multi-objective model. 


optimize Z^ 

1 — ■ 1 j • • * | 

s 

S • *fc 9 

AX < 

b 


X > 

0 


Adding aspiration levels and deviations variables to each 

objective and weighting 

each resultant goal, 

s 

we obtain, 

Minimize a = 

Z (w-a - + r- p. ) 
i=l 11 11 


subject to 


Z^x) + - « ± = 

AX < b 

Z° , i = 1, . 


x > oc ± , > 0, 

I 25 1 f » * « | s 


where. 


a i = Weighting factor for the positive deviation 
of goal i. 

(3^ = Weighting factor for the negative deviation of 

goal i. 

Z? = Aspiration level for objective 1 . 

Interval^ Go al_ Prog ramming j 

In this model the decision maker specifies a range 
of aspiration levels for each of the goals instead of one 
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aspiration level* It obtains ar optimum solution by minimiz— 
ing the weighted sum of deviations from the set of ranges. 

The mathematical formulation of this model is, 

s 


Minimize a 
subject to 


= T 
i-1 


(w i, 2 “i, 2 + r i, 1 P i, 1> 


2 


i,i ' ^ 1,1 “ a i,: 


^i,2 J '" 2 " a i 2 


L 1 all s 

U i all s 


where. 


Z i,l = Z i,2 

= expression for objective k. 

^i,l» p i,2 

= negative deviatL ons . 

a i,l* a i, 2 

= positive deviations. 

r i,l 

= weighting factor for the negative 

deviation for goal Z. .. , 

1,1 

C\J 

d‘ H 

= weighting factor for the positive 

deviation for goal Z 

i,i 


Fuzzy Ljnear_ Programming ° 

A fairly recent attempt at modeling and solving the 
multiple-objective problem is that known as fuzzy programming. 
The approach is similar, in many respects, to the weighted 
linear goal programming method previously discussed, differing 
primarily m the manner in which the importance of the goals 
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are considered. This method minim' res the worst under achieve- 
ment of any goal. 

A major advantage of fuzzy linear programming is that 
it may be transformed into a conventional linear programming 
model. The main disadvantage of this method is, the under- 
achievement of just one goal can have a major impact on the 
solution, since it attempts to minimize the maximum under- 
achievement. 



CHAPTER III 


NETWORK GOAL PROGRAMING ALGORITHMS 

3.1 ^troduction » 

In this chapter, we shall consider two special classes 
of linear goal programming problems, i.e. weighted goal network 
flow problem and interval goal network flow problem. 

We generally confront problems that require two objectives 
to be considered simultaneously. Often, these objectives are 
conflicting m nature. These types of problems may arise in 
road networks, co m munication networks and pipe line networks. 

The typical applications of bi- criteria network flow problems 
are mentioned in Section 1.1. 

One of the methods of solving bi-criteria problems 
is to obtain the efficient or non-dommated solutions set. 

As the efficient solution set is considerably large, it. becomes 
practically difficult for thest decision maker to choose the 
solution he would prefer. 

Goal programming is a technique which takes the deci- 
sion maker' s preference into consideration and provides him 
the preferred solution. He expresses his preference by 
specifying targets to the objectives and the weights to the 
deviations from the targets. 
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In weighted goal programrn lg, the decision maker speci- 
fies a target for each of the objective and weights for the 
deviations. It provides him an opcimum solution which is as 
close as possible to the specified targets, by minimizing 
the sum of weighted deviations. 

In interval goal programming the decision maker speci- 
fies a range or interval for each of the objectives and weights 
for each of the deviations from lower and upper bounds of the 
interval. It obtains an optimum solution lying in the ranges 
specified if one exists or finds a solution such that the sum 
of weighted deviations are minimum. 

This chapter has been divided into 11 sections. In 
Sec. 2 we present the notations used to represent the feasible 
region of the bi-criteria network flow problem. The tracing 
of various trade-off curves is presented in Sec. 3. The 
mathematical formulation and development of the WGNF algorithm 
are presented m Sec. 4 and 5. In Sec. 6 and Sec. 7 we present 
the statement of the "WGNF algorithm and a numerical example. 

The mathematical formulation and development of IGNF algorithm 
are given m Sec. 8 and Sec. 9. 4 numerical example of IGNF 

problem is given in Sec. 10. Finally the computational perfor- 
mance of both the algorithms is reported in Sec. 11. 
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3.2 Notations^ 

The various regions as shown in Fig. 3.1 are defined 
as below; 


E ; 

{ (x,y) 

p 

« 

x > 0 and y > 0 }. 

S $ 

{ (x,y) 

a 

x = 2 c. x. , y = 2 

(i,o)eA (i,o)eA 

'1= 

{(x,y) 

0 

0 < x < f-^, 0 s y g 2 and (x,y) ) 


{ (x,y) 

V 

0 < x < £y g 2 < y < « and (x,y) ) 

D 

{ (x,y) 

C 

£3 \ x < °° 9 84 < y i a and (x,y) j 


{(x,y) 

o 

f l i x < », 0 < y < g^ and . x,y) j 


{ (x,y) 

£ 

Si x < f-J 

s 2 ; 

{ (x,y ) 

8 

Si x < f,' 

3 ' 


x 


„}• 


Sy, * 


a l* 

h' 


V 


w. 


{ (x,y) e S ; x > 

{(x,y) e S ; x > f^} 
positive deviation from goal i, 
negative deviation from goal i. 
weight assigned to positive deviation a 1 * 
weight assigned to negative deviation . 
r, > 0. 


'i» *1 

3.3 Tracings the^ Trade-Off Curves 


The optimum solutions for VGNF problem and IGNF problem 
lie either on B^, B^ and B^ curves or in the feasible 
region S as shown in Fig. 3 . 1 . Since our algorithm searches 
for an optimum solution by traversing along U II B^ U B^ 



Fig. 3 


W) 


I 

I 

\ 

*2 ! 




1 feasible region of b»cr it^ria network 
flow problem 
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and through the feasible region S it is necessary to trace 
the B^, B£, and B^ curves. Once these curves are obtained, 
it can be shown that they to-getner form a closed convex fea- 
sible region S, 

Tracing B^ Cu rve s 

The detailed procedure for obtaining the B^ curve was 
presented m Sec. 2,4. The procedure for tracing B^ and 

B^ curves is same as that for tracing B^, except that the 
optimality conditions, the critical set, the active set, the 
passive set and the criteria for selecting the price ratio 
differ, as the problems considered to obtain each of these 
curves are different. 

Tracing B 0 Curves 

B^ curve is a piecewise linear concave function between 
cost and budget. This curve is obtained by parametrically 
solving the following problem. 


Max. Z = I c . x. 

(i,j)eA ^ 

subject to 


r 

\ 


-V, if i=l 


2 x . - 2 x. - « « V, if i=n, V ieN 

(j,i)el(i) 3 (i, j)eO(i) x3 / 


L- 0, otherwise 


0 <. x. < b . 
- lj - ij 


(i, j)eA 


V (i, j) e A 
d i0 x ij 4 D 
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The optimality conditions are 5 
(±) “pq-^pq - ° 

(ii) -c. + XL. S’. . 0 V (i, J ) e L U u 

-L J _L J 

(ili) UX D- 2 d . x . _ ) = 0 

(i,o)eA 

The non-basic arcs are classified as s 

(i) Active arcs S' $ {(i,o) e L U Us. c. < 0 and d. _ < 0} 

IQ J- J 

(ii) Critical arcs S":{(i,j) e L U Us c 0 and d. < 0} 

(iii) Passive arcs S ,M ;{(i,;j) e L U Us d >0} 

f 0 

Initially, £ c. x. is maximized and let 
U,0)eA ^ ^ 

(C , D) be the solution, A critical arc (p.q) e S " is 

IIlcXA. * 

selected and entered into the basis and dual-simplex iteration 
is performed. This is continued until S is empty. 

Next, an active arc (p,q) that that '(-A = min. , (^_) 

pa - (i.a)es' 

is selected and entered into the basis and dual-simplex iteration 
is performed. Each dual-simplex iteration results m an extreme 
solution. This is repeated until S' = {0}, As ti is gradually 

increasing and ultimately becomes zero along with gradual 
decrease m C and D, the curve is a piecewise concave linear 
function. 

Tracing^ Curve - 

is a piecewise linear concave function between cost and 
budget. This curve is obtained by parametrically solving the 
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following problem. 


max. Z s= E c. x. . 
(i,j)eA 


subject to 

2 x 

(j, i)el(i) 

0 < x. < b 
“ 10 " 

2 d, 
(ifj)eA 


O 1 


xj 



f-V, if 1=1 

2 x n = ) V, if i=n, ¥ ieN 
(i,j)eO(i) /^o, otherwlse 

¥ (i,j) e A 

> D 


The optimality conditions are* 

(i) o pq+ X<d pq = 0 

(ii) - c l0 -/td iG >0 ¥ (i,j) s A 

(iii) -£C( 2 d. . x. - D) =0 

(i,j)eA 

The non basic-arcs are classified asi 

(i) Active arcs S' i {(i,j) e L U U; c\ < and d. > 0} 

i J XJ 

(ii) Critical arcs S "a {(i,j) e L U Ui c. ^0 and d - . > 0} 

X J "* X J 

(iii) Passive arcs s'"s {(i,j) e L U U; d. < 0} 

X J 

Initially, 2 c. x. _ is maximized and let (C , D) be 
(i, j )eA ^ ^ max 

the solution. Next, the critical arc set is emptied by performing 

dual-simplex iterations. An active arc (p,q) such that 

= min. ( ^40 is selected and entered into the basis. 
pq (i,j)eS 

The dual simplex iteration is performed. This is repeated until 


\ 
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S = {0}. Since the value of /t Is gradually increasing and 

C is gradually decreasing with the increase m D, the curve 
is a piecewise concave function. 

Trac ing Curvey 

curve is a piecewise linear convex function between 
cost and budget. This curve is obtain edh by parametrically 
solving the following problems 


Min. Z = 


1 c x. 
(i,o)sA 


s.t . 


C-V, if i=l 

1 x . - I x.. =; V, if i=n, V i e N 

(j»i) £l (i) 0 (i,j)eO(i) ^ 


0 < x. 

•’ 10 


< b. . 

10 


O f otherwise 
¥ (i,o) e A 


Z d. x. . > D 

(i,0)eA ^ ^ 

The optimality conditions are* 


( i ) c -lid =0 

pq pq 

(ii) c ± j - Atd.^ } 0 


V (i,o) £ L U u 


(lii) -46 ( I d. . x. - D) = 0 

(i,d)eA ^ *3 

The non-basic arcs are classified as? 

(i) Active arcs S' s i (i, j) £ L U U s c. > 0 and d. > 0} 

X J 1 J 

(ii) Critical arcs S f h -(i»o) e L U U i c. . <, 0 and d. > 0} 
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(iii) Passive arcs S ,M ; {(i,j) e L U U s d. 1.0} 

X J 

Initially 2 c x. is minimized and let (C^^, D) 
(i,o)eA 13 1 0 

be the solution. Next, the critical arc set is emptied by 

performing dual-simplex iterations. An active arc (p,q) such 

that lb = mm. (-^O is selected and entered into the 
pq (i,o)eS’ 

basis. The dual simplex iteration is performed. This is repeated 
until S' = {p}. A dual simplex iteration results in an extreme 

solution. Since the value of IL is gradually increasing and both 
C and D are increasing, the curve is a piecewise convex 
function. 

The curves and B^ together constitute a convex curve 
as the starting solution is same and the cost is gradually 
increasing. The curves B^ and B^ together constitute a concave 
curve as the starting solution is the same and the cost is 
gradually decreasing. It follows from the above discussion 
that the terminating solution of B-^ is the starting solution 
for I$2 and the terminating solution of B^ is the starting solu- 
tion for B^. Similarly the terminating solution of is the 
starting solution for B^ and the terminating solution of is 
the starting solution for B^. Thus we obtain a closed region 

enclosed by the curves U B 2 U B 3 U B 4* Since B i u B 4 is a 
convex curve and B£ U B^ is a concave curve, the closed feasible 

region xs & closed convex region. 
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3*4 


Mathematical Formulation of the Weighted. Goal Network 
^owWoFlemV * * 


In this section, we give the mathematical formulation 
of the weighted goal network flow problem. 

In the network G = (N, A), associated with each arc 
(i,j) e A are two numbers c. and d . The capacity of each 

1 J 1 J 

arc (i,j) e A is b , The amount that is to be shipped from 

1 J 

source 1 to sink n is V. 


Minimize Z = w-^a^ + r-^P-^ + w 2 a 2 + r 2$2 ( 3 • 1 ) 

subject to the following constraints i 


Flow_ Cons ervat ion Constraints* 

These constraints essentially represent the factt that 
flow of the commodity is conserved at all nodes, except at 
source and sink. 


I x . 
(j,i)el(i) J1 


E x . 
(i,j)eO(i) 10 


i 

_ ) 


L 


-V, if i=l 


V. 


0. 


if i=n, V i e N (3.2) 
otherwise 


Capacity Constraints 5 

« m -■j Uu.nc.-aa w L.'mii rm-rmk -*rr* m j* tax . »- a 

These constraints essentially represent the fact that 
flow over an arc (i,j) cannot exceed its capacity 


0 < 


X i0 


< b 


10 ’ 


V (i,j) e A 


( 3 . 3 ) 
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Goal, Constraints z 

An objective in conjunction with an aspiration level 
and deviations is known as goal constraint . 


, 2 c. x. . + p-, - a n = a 

(i,j)eA ij ij 1 1 1 


, 2 N „ d x •+ P p - a ? = D 

(i,j)eA ij 2 2 1 


(3.4) 


(3.5) 


3 . 5 Devel opment of ^ the_Algpr ithm_i 

We shall now discuss, the various steps involved in 
finding an optimum solution for a specified (C^, D,, w-p r^, 
w 2* r 2^‘ 4et us c °nsider ( C ]_» D ^) £ R 2 in ord - er "trace all the 
possible steps of the WGNF algorithm. 

Initially, the minimum cost flow problem with 2 c._x., 

(i,j)eA 10 ^ 

as an objective function is solved. Let (C ^ n , D') be the minimum 
cost and budget obtained and it is represented by I in Fig, 3.3. 

The deviations axes from aspiration levels (Cp D^) are shown 
in Fig. 3.2. If the algorithm traverses in quadrant I then it 

minimizes and a 2 and the other two deviations are zero. 

Similarly and ^ are minimized in quadrant II, p 2 and 

are minimized in quadrant III and « 2 and |S^ are minimized in 

quadrant IV. Therefore, the algorithm always minimizes only 
two deviations and the other two being at zero. The above 
statement is supported by theorems 1 to 5 given below. 




Fig. 3-2 Aspiration tevets and devaiations for 
W6NF problem 


■ II i I 
i 1 (C V 0, ) 



Fig. 3 3 Paths traced by WGNF regorithm to obtain 
optimum solutions 
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THEOREM IS = 0, V i = 1,2 

Proof* We will prove this theorem by contradiction. 
Suppose (X*, a^, p£) be an optimum solution and ct£, p^ > 0, 

Let C* = I c. xf 
(i,j)eX 


The value of the objective function is 


\ = w^a* + r x pj (3.6) 

The goal constraint is 

C* + Pi - «i = c ! (3.7) 

Let A = min. (p^, a-f ) 

Consider the solution (X* , p^, a^), 

let. Pi = P* - A (3,8) 

a l “ a l “ A (3.9) 


substituting (3,8) and (3.9) in (3.6) and (3.7), we get, 

= w^(a^ + a) + r^(p^ + a) 

= w^a-^ + + a(w^ + r^) (3.10) 

goal constraint is 

C * + p l “ a l = C 1 (3.11) 

Hence the goal constraint is satisfied. 


Since A, w-^, r-^ > 0 it follows that > Z^. This 
contradicts that (X , p^, a-^) is an optimum solution. Therefore 
a = 0, Hence the theorem follows. 
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Similarly it can be proved lor a 2 P 2 “ Therefore 
atleast two ol a-^, a 2 and P 2 are zero. 


QED 

THEOREM 21 If (C^ 3^) e i^, then (C*, D*) e B ][ and 
D* and C* > c^. 

Proof; Consider (C D ') e S 2 > There exists (O', D) e S- L , 
such that D < D Since D < D' , it follows that a 2 < a 2 resulting 
in Z < Z ' where Z is the value of objective function (C’ } D), Hence 
a better solution exists in S 1 for every solution in S 2 . 

Consider a feasible solution (C, D’) e S^. A feasible 
solution (C, D 1 ) e exists such that C < C', Since C < C’ it 
follows that < a^« ft results in Z < Z' where Z and Z 1 are 
objective functions of (C, D ') and (C 1 , D') respectively. There- 
fore for any solution in S, a better solution exists on B^. 

Let (C 1 , D') e B^ and D' < D^. There exists (C, D-^) e B^ 
and C < C'. Therefore < a-^ resulting m Z < Z 

Let (C ", D ") e B^ and C " < C^, There exists (C-^, D) and 
C-, > C " , Therefore a 2 < a 2 " and resulting Z < Z 

Hence an optimum solution is (C*, D*) e and D* £ D, and 
C* > C-^. It follows from theorem 1 and above that if (C^,D^) e R^ 
then the optimum solution is (X*, a^, a 2 ) and (3-^, (3^ = 

Similarly we can prove the following theorems, 

GKD 
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If ^ C l» V e R 2 then D *) e B 2 and D *^ D l 

and C * < C, . 

The optimum solution is (X* a 2 , Pp) and a-j_, P 2 = °« 

If ^ C l» V £ R 3» ther " ( C *» D *) £ b 3 D * < D ! 

and C* < 

The optimum solution is (X*, ( 3 *, p *) and a^» a 2 = °* 

If (C lt D 1 ) e R 4 , then (C*, D*) s B 4 and D* < ^ 
and C* } C^. 

The optimum solution is (X , a-j_) ana a 2 , P-j_ = 

OBD 

If < D' then the algorithm moves along the 
curve in the direction of J as shown in Fig. 3.3 ► Otherwise 
it traverses along B 4 « The algorithm evaluates the objective 
function values of extreme- solutions and two non-extreme 
solutions, namely (C, D^) e B-, U U Bj U and (C^, D) a B^ U 
B 2 U U B 4 , because one of these solutions is optimum if 
(C^, D^) ft S. ¥e ascertain the above statement in theorems 6 
to 9 given below. 

THEOREM 6; If (C-^ D^) / S, then one of the extreme solutions 
of B 2 U B 2 U U B 4 is optimum. 

Proof i Consider two adjacent extreme solutions (C 1 , D ') and 
(C ", D ") . Any non- extreme solution (C, U) lying between 
(C D ') and (C ", D ") can be expressed as convex combination of 


THEOREM 3 s 


THEOREM 4s 


THEOREM 5i 
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(C , D ) and. (C , D ), Therefore Z 1 or z" is atleast as good 
as any Z, depending on Z' Z" or Z" Z'. Hence the 
theorem follows. 

QED 

THEOREM 7* If (C.^, £ lies between two adjacent extreme 

solutions (C\ D ') e and (C ", D") s such 
that C-^ > C ", D 1 > and satisfying Z' > Z where 
Z', Z and Z" are objective function values of 
(C', D'), (C^, D) and (C ", D ") respectively, then 
< C 1- D) is optimum. 

Proof; Z " = ^2 a 2 + r iPi (3.12) 

Z — w 

Since < a^» it follows that Z < Z". (C^, D) is 
atleast as good as any other solution lying between (C-^, D) and 
(C 1 , D ') as they can be expressed as a convex combination of 
(C.^, D) and (O’, D') and Z ' ^ 2. Similarly (C^, D) is better 
than any other solution lying between (C^, D) and (C ", D ") because 
they can be expressed as a convex combination of (C-^, D) and 
(C ", D ") and Z" > Z. According to theorem 2, for any solution 
m S, a better solution exists on . Hence (C^, D) is atleast 
as good as any neighbouring solution. Therefore (Ci , D) is a 
local optima. In linear programming, local optima is also a 
global optima. Hence the theorem follows. 
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Similarly this type of result can be proved for each of 
(Ci» D) 6 Dp) £ and (C, D^) e lying between two 

adjacent extreme solutions, (C D' ) and (c", D ") satisfying 
2 ' >. Z then they are optimum, 

QED 

THEOREM 8i If (C, D^) s B^ lies between two adjacent extreme 
solutions (C 1 , D') e and (C ", D ") e such 
that C > C ", C " > and satisfying 2 < z" where 
Z, Z', Z " are objective function values of (C, D^), 
(C, D’) and (C ", D") respectively, then (C, D^) is 
optimum. 

Proofs Z' = w p a i + r 2^2 
Z = w-^a^ 

Since a p < cep > it follows that Z <; Z\ (C, D^) is a 
better than any solution between (C, D') and (C, D-^) and 
between (C, D^) and (C ", D") as they can be expressed as 
convex combination of these solutions and Z < Z ' and Z < Z ", 
According to theorem 2, for any solution in 3, a better solution 
exists on 3^« Hence (C, D^) is atleast as good as any neigh- 
bouring solution. Therefore (C, D^) is a local optima. In 
linear programming, local optima is also a global optima. Hence 
the theorem follows. 

Similarly this type of result can be proved for each of 
(C, D-^) e B^, (C lf D) e B £ and (C^, D) e B^ lying between two 
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adjacent extreme solutions, (C', D') and (C", D") satisfying 
2 < 2 " then they are optimum. 

QED 

THEOREM 9 * If (C, D^) e and (C^, D) a lie between two 
adjacent extreme solutions (C, D ') e B 1 and 
(C ", D ") e such that > C ", then 
U) (0 1( D) is optimum if Z <, Z and 
(ii) (C, D^) is optimum if Z ^ Z 
where Z and 2" are the objective function values 
of (C lf D) and (C, D^) . 

Proofs ¥e shall prove case (i) . 

Let Z' and Z " be objective function values of (C 1 , D ') 
and (C", D") 

Z ' = w 2 ap r^^ 

2 = w 2 a-2 

Since a 2 < a 2 , it follows that Z < Z". (C-^, 0) is 

atleast as good as any solution lying between (C^, D) and (.C, D^) 
because they can be expressed as convex combination of (C-^, D) 
and (C, D^) and Z <_ Z, Similarly (C^, D) is better than any 
solution lying between (C^, D) and (C ", D ') as they can be 
expressed as convex combination of (C^, D) and (C ", D") and 
Z < Z". According to theorem 2, for any solution in S, there 
exists a better solution on Hence (C^, D) is atleast as 
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good as any neighbouring solution. Therefore (C^D ) is a 
local optima. In linear programming, local optima is also a 
global optima. Hence the theorem follows. 

Case (ii) is proved as follows i 
Z' = 
t = 

Since a ]_ < » it follows that Z < z' using theorem 2, 

Z < Z T and Z £ Z, it can be proved that (C, D^) is a local 
optima, which is also a global optima in linear programming. 

Hence (C, D 1 ) is an optimum solution. 

Similar results can be proved for the cases when 
(C, D 1 ) e B 2 U B 3 U and (C^ D) e B 2 U B 3 U B^ lie between 

two adoacent extreme solutions then one of them is optimum 
depending on Z > Z or 2 < Z, 

QED 

We can infer from theorems (7) to (9) that the non- 
extreme solutions (C, D 1 ) and (C ]L , D) behave as extreme solutions 
and hence the algorithm evaluates the objective^ function values 
of these non-extreme solutions. 

If Z corresponding to (C, D) the current solution is 
greater than or equal to Z " of (C ", D") the previous extreme 
solution, then (C", D") is optimum and the algorithm stops. This 
is proved in the following theorem. 
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THEOREM 10! If (C^ e R 1 and (C, D ') e (tf, ff) e B 1 
and (C ", D") e B^ Be three adjacent extreme 
solutions satisfying Z" > Z and Z'>^ Z then 
(c; d) is an optimum solution. 

Proofs (^, E>) is atleast as good as any other solution lying 
between (£T, IT) and (C 1 , D ') as they can be expressed as convex 
combination of these two solutions and z’>^ Z. Similarly 
(C, D) is atleast as good as any other solution lying between 
(S, CO and (C*', D") because they can be expressed as convex 
combination of these two solutions and Z " > Z. According to 
theorem 2, if (C^, D^) e then a better solution exists on 
Bj than any solution in S. Therefore (C, D) is a local optima, 
as it is better than any neighbouring solution. In linear 
programming local optima is also global optima and hence 
(C, EO is optimum. Similarly, this type of result can be 
proved for the cases, (C^, D^) e R 2 U U R^. 

OED 

If <_ D n , then the algorithm moves along B^ U 
and obtains an optimum solution. Otherwise the algorithm 
halts at (C, D-^) e B^ represented by J m Fig. 3.3. If (C-jy ^ 
e R^ then (C, D^) is an optimum solution which follows from the 
theorem 8. Otherwise it moves vertically up along <£L in the 
direction of C^. 
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We will discuss in detail how the algorithm increases 
the total cost C, while maintaining the total budget D constant. 
This is achieved by simultaneously augumenting the flow in two 
cycles. 


An arc (i,j) is a forward arc if the flow is sent in 

the direction of the arc, and a backward arc if the flow is 

sent against the direction of the arc. We will refer to this 

as orientation of the arc. The basis structure of current 

solution is (T U {(p,q)}, L, Uj and its basis consists of 

exactly one cycle formed by the arc (p,q). This cycle is 

denoted by W . Now, define a set 
PQ. 

S° ={(i,c) eLUU: dpq/to < ° 3214 



d 

- 


10 



> 0 ) }. 


The above conditions are derived as follows. 

Criteria for Selecting, a Non r Basic_Ar,c for Z Cycle, flow 
lugumenjatip >nV~" 

We will derive the conditions that the non~basic arc 
(r, s) e L U U must satisfy, such that by simultaneously augument- 
ing flow in cycles W pq and \l rs , the total cost increases and the 
total budget remains constant . 

The basis structure is (T U {(p,q)}> L, U) . 
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Le-t A pq be tlie flow bo be circulated, in cycle W . 

pq 

e "^ A rs bo he circulated in cycle W 

J rs 

For the total budget D to remain constant, 

A pq kpq + A rs d rs = 0 (3.13) 

In order to increase the total cost C, 


A c 

pq pq 


A c 
rs rs 


> 0 


Using (3.13), 



d 

( : EiL ) 


= K 


rs 


(3.14) 


(3.15) 


Since a 


pq 


A 


rs 


we get. 


0, hence m order to satisfy (3.13), 


JE>£. 

*rs 


< o 


Dividing (3.14) by a 


pq 


pq 


d. o 

■Ea x c 

**XM 

d rs 


rs 


and substituting (3.15), we get, 

> 0 


( 3 . 16 ) 


(3.17) 


Hence if an arc (r, s) satisfies (3.16) and (3.17) and augumenting 
flow a^ and A^^Kin cycles VT and ¥^ s , the total cost increases 
and total budget D remains constant. 

If S° = {^} then reverse the orientation of the arc (p,q) 
and set c pq =~ c pq , d^ = - d pq . Then, we again find S°. If S° 
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is not empty , select an arc (r,s) e S° and. enter it into the 
basis* This results in the formation of two cycles W and 
W rs* ^ orientation of a cycle W. (i # j) e L U U is same as 
that of the arc (i-, J ) • Let and be the sets of ares, 
along and opposite to the orientation of the cycle Vh , res- 
pectively. Let A be the amount of flow to be augmented in 
V/ and K.a be the flow to be augmented in cycle W where 

Jr St Jr M. X O 


K 


- Spq/cfpg* Next, we determine A as follows. 

Consider an arc (i,o) e (W fl W _)* The flow in arc 

Ti rl X S 


pq 


(i,j) after augmenting flow in W and W is, 

pq rs 

x. • = x. . + a L + L K A 

10 10 pq 1 2 pq 


where , 


id 


" x ij + A pq ^1 + I 2 


=s existing flow m (i, j ) 

= +1, if (i»d) e K 


pq 


- 1 , if (i,o) e W 


pq 


= + 1, if (i,d) e % 
= - 1, if (i»0 ) & 1 


rs 


-rs 


A (1,+IJC) = resultant augmented flow in arc (i,d)* 

pq 1 2 

In order to satisfy the capacity constraints, 


0 < x. . + a 

- io pq 


(I 1+ I 2 K) 





50 


If (II + I 2 K) > 0, then 


a < 

pq 


(b . ~ x ) 

(i-j+I^) 


If (1^ -! I^K) < 0, then 


A < 

pq - 


fij, 

(Il+I 2 K) 


In this manner A pq is calculated for all arcs 
(i,j) e (V pq U U rs ). Therefore, 


A 

pq 



X. 

min {- if f I -+T K'i <• n 

; 1( o) E (W pq UW rs ) (I 1+ X 2 K) 1 +I 2 K) 

(b u ” x n> 

~~* L - if In+IpK > 0} 

(i 1+ i 2 k) 1 2 ' 

, (C.-C) 

min. (a , . , ) 

pq (c + Kc ) ^ 
pq rs y 


Once A pq is obtained, the flows are updated in ¥ pq and ¥ rg cycles. 
One of the arcs (u,v) e (W pq U W rg ) leaves at one of its bounds 
respectively. The dual variables, c 1 s, d\ * s, V (i, j ) e L U U. 
C, D and Z are updated respectively. 


Arc (p,q) for the next iteration is obtained as follows? 


(a) If (u,v) = (i,o) £ (¥ pq ) , set (p,q) = (r,s). 

(b) If (u, v) = (i f j) e ({¥ pq 0 ¥ rs } U Wj, (p,q) does not 
change . 
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If "the arc (p,q) is at one of its bounds, reverse the 

orientation of (p,q) and set c = -c. '& n =-d . By select- 

PSl PP Ph 

ing an arc (r, s) e S at each iteration and augmenting flows 
ln ^rs* a lS°rithm moves along J,L shown in 

Fig. 3# 3. If £ 3, then the algorithm obtains this 

solution which is an optimum solution and hence stops. If S° 
is empty then the current solution, 

(C, D) e if °Pa > 0 

0r » 

(C, D) e B, if < 0 

D d 

pq 

according to the theorem 11. L in Fig. 3.3 represents this 
solution. 

THE0REI1 11; If S° = ^ 0 } then the current solution (C,D) e 

WV > ° and (C,D) £ ®3 i£ W'Vq “ °‘ 

Proof; The necessary and sufficient conditions for (C, D) e B 2 
are, the >non~basic arcs ( 1 , 3 ) e L U U whose > 0 and d ^ 0, 
do not exist. 

The necessary and sufficient conditions for (C, D; e 
are, the non-basic arcs (i,j) £ L U U whose b in > 0 and d i;) >, 0, 
do not exist. 

Suppose, > 0 


4 5’* * Aft 

ft I a ,**.* , A riBt I 


mm mm m? 


8747 Pi 

jfctk mu -Mfe* 
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Arcs of the type c_^ > 0 and d^ < 0 satisfy ( 3 . 16 ) and 
(3.17) if Cpq > 0 and d^ > 0 , Arcs of the type c^ > 0 and 
d. = 0 satisfy ( 3 . 16 ) and (3.17), if c x 0 and d y 0 . 

pq pq 

Since S = l 0 ; , the non-bo sic arcs ( 1 , 3 ) of the typ>e 
>c i 3 0 and ^io • < . 0 are absent. Therefore (C, D) e B 2 « 

Suppose, c /d <.0 
pq pq * 

d 

Arcs of the type c. > 0 and d. > 0 satisfy < 0 and 

■ d 13 ' 13 d. 

c pq - ■ Pa X b ij > 0 lf °pq. > 0 ^ V ' °- S ° - ' 0 } 

±0 

the non-basic arcs ( 1 , 3 ) whose c > 0 and d\ ^ >, 0, 3f 
v(i, 3 )eLUU are absent. Hence (C, D) e B^. 

QEO 

Now, the algorithm moves towards M along the curve B^. 
According to theorem 2, if (C^,D^) e the (C , D ) e and 
D* } D-^,and C* < C^, it is justified m moving along L,M for 
finding an optimum solution. 

Similarly if (C^, D^j e R^, the algorithm traverses along 
I, F, H, Y and finds an optimum solution. If (C^, D^) e 
and D 1 > D ffiax , then the algorithm moves along and curves, 
and obtains an optimum solution. 
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3.6 Statement^ of Algorithm". 

A. formal statement of the algorithm organised m a 
manner suitable for computer implementation is given belov/5 

Step. 1* Set count to zero. Solve the minimum cost flow problem 

with I. c- x as the objective function value. 
(x.l.eA ^ ^ 

Let (B,L,U) be the optimum basis structure obtained 
and X be the solution. Compute, 


C as Z C X 

(i,j)eh 10 1J ' 


and D = 


Z d . x , 
(i,j)ea 13 13 


(C, D; is the minimum cost and budget respectively. 

Treating B as the basic tree, define the dual variables 
c d. 

ix 's and ix 's as follows. 

J J 

U 1 C = °» = 11 1° + c ij » ¥ (i.d) e E 

•j, ^ = 0, ji d - ix . d + d. ¥ (i,l) e B 

X 0 X xj 

Define c and d. for all non-basic arcs (i,o) a L U U 

10 1 J 

as follows s 


c 


xo 





-I- c. 


xo 




¥ ( i, o ) e L 

¥ (i,l) £ U 

¥ (i,j) e L 
¥ (i,o) e U 
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Comput 

e y 

^ (C - 

i — i 

o 

* 



a 2 

= (D - 

- V 



Pi 

- (C 1 

- c) 



p 2 

= (D, 

" 0) 



\ a r 

if 

> o, 

v i = 1,2 

a, = 




i 

‘ o, 

if a 

s 0, 


p i = 

K r 

if p ± 

> 0, 

¥1=1,2 

l°> 

if q 

< o, 


Z = 

w 1 a 1 + 

w 2 a 2 + 

r l p l 

+ r 2 p 2 

If D, 

< D then flag 

= “1, 

else flag = 1. 


Set Z = Z, X = X, c' = C, D* = D and go to Step 2. 

Step. 22 Let S' = £ L U II: c. . <0 and flag *d. > 0} 

J _L J 

If S ' - 0 then go to Step 3* Otherwise select an 

arc (p,q) £ S' and enter it into the basis. Perform 
pivot iteration, up date x. ' s. 

-L J 

If flag * D > flag * D-, , then circulate A flow against 

1 P-D 

the orientation of 'I where a = , and go to Step 7. 

J/H. (J 

, , Pq 

If Z > Z then X is an optimum solution, Stop, otherwise, 

update (B,L,U), a c, s, ft ^'s, set Z ' = Z, X' = X, C 1 = C 

J J 

and D' = D. Repeat this step. 
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Step. 2 s Let S' = 
Define, 


e L U U; c. > 0 and flag * d. > 0} , 

- 1 - J X J 


% - ■^ i * £la s- 

id 

If S = \f>jj then set flag = - flag and go to Step 4. 

Otherwise select an arc (p,q), M = mm. (44), and 

Pq (i,o)eS’ 13 

enter it into the basis. Perform the pivot iteration, 
update x ’s. 

1 J 


If C > > C ' then circulate A flow against the 

orientation of W where a = (C-Ch )/c . Set Z' 

PI 1 ' pq 

X' = X, C ' = C and D' = D, Augment a flow m T J 

pq 

A = (0- Cp/O^. 


= Z, 

where 


If flag * D > flag * D^, then circulate A flow against 
the orientation of W where a = (D-D.. )/d , then check 

Jr tL “*“ ir ri 

if Z Z' then X' is an optimum solution, stop. 
Otherwise, go to Step 7. 


If L y Z ' then X' is an optimum solution, stop. 

c d, f 

Otherwise update (B,L,U), n. 's, tc s, set Z ' = Z, 

J J 

X' = X, C' = C, and D ' = D. Repeat this step. 

S£ep. 4s Let S' = [(i,j) e L y U; c^ > 0 and = 0 ) . 

If S ' = 0 then go to Step 5. Otherwise select an arc 
(p,q) e S' and enter it into the basis. Perform the 
pivot iteration, update x. 's. 
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I f ° > C i > c ' then circulate a flow against the 
oiientation of ¥ where A = (C— C-,)/c' » X is an 

ir H -1 pC[ 

optimum solution, stop. 

If Z Z then X 1 is an optimum solution, stop. 

Otherwise, update (B,L,U), % 0 - s , n d < s , Set Z' = Z, 

v , / f J 3 
X = X, C = C and D = D, Repeat this step. 

Step. Let S' * i(i,j) e L U U; c > 0 and flag * d. > 0] . 
Define, ^ •flag/'d^. 

If S’ = 10} then X' is an optimum solution, stop. 

Otherwise select an arc (p,q) - max ) 

pq (i,a)eS’ 

and enter it into the basis, go to Step 6. 

Step. 6; Perform the pivot iteration, update x ' s. 

J 

If C > > C' then circulate a flow against the 

orientation of V where A = (C-C^/c^. Then if 
Z y Z’, X' is an optimum solution, stop. Otherwise 
X is an optimum solution, stop. 

If Z > Z' then X' is an optimum solution, stop. 

Otherwise update (E,L,U), tc c 's, x d, s. Set Z ’ = Z, 

J J 

X' = X, C = C and D' = D. Repeat this step. 

Step Jr- If count = 2 then set Z = Z, X = X, C = C, D = D 
and go to Step 8. 

If C /. C lt then X is an optimum solution, stop. 
Otherwise the basis structure is U { ^p>q.^ j jL,U), 
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Define S° = t(i, D ) e L U Us d /d. < 0 and 

PQ. 


Cc pq ~ ( W d i3 ) * C id > 0) )- If s ° is 


empty then set count = count + 1, c = -c 

5 PI pq* 

^pq = status = -status and go to Step 7, 

Otherwise, select an arc (r,s) e S° and enter it 
into the basis (T U {p,q}). 


Define, 


L . 

ij 


i are formed i. 

e - l pq’ '“rs- fle£Lne 

/d’ , set 

( rs ’ 

T. 

11 

= 0, V (i, 3 ) e A. 

i ’«“■ 

if 

(i,j) e W 

pq_ 

/ 7 -1, 

if 

(i - 3) e *pq 

( 7 +K, 

J 10 

if 

(1,3) £ V rs 

/ T. -K, 
[ U ’ 

if 

(i,3) a y TS 



X . 


A = 


H 


min. (- if T < 0, 

(i.oWV'rP Id 13 

) t if t > o). 

±3 ° 


pq 


(c, -c) 

min. (a, J 

(c +K* c„_) 

pq rs 


Update A p(J in » n0 and K 4 ^ in W rg . 


pq 


pq 
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Define S° = {(i,o) e L U Us ^ /Si^ < 0 and 

(c pq > 0) - If S ° ±S 

empty then set count = count + 1, c = -c . 

pq pq 

= -d , status = -status and go to Step 7. 
pq pq ° ^ 

Otherwise, select an arc (r,s) e S° and enter it 
into the basis (T U {p,q}). 


Two cycles are formed i.e. W , ~<‘ rs * Define, 
K = -<CAi set T ia = °> V (i,o) e A. 


pq rs 


Define, 


T. = J 
U 


( qyi, if (i,a) S w 


pq 


(V 1 ' if (l ^ )e¥ pq 


T. . = 
10 


\ T +K, i f (i» 0 ) e ¥ 


i 

/ T 
l ^ 


rs 


-K, if (i,j) s 17 


TS 


A = 


A 


pq 


mm. (- if T < 0, 

(l ^ )£(W pq U 'r 5 ; « 

(hi “«■ ), if I > 0). 

id 3 

(C, -c, 

= min. (a, - ■ ) 

(c -:-K * c’ ) 
v pq rs 7 


Update a in ¥ and K a m ¥ « 

* pq pq pq rs* 
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Update (B,L,U), % C 's, % d ' s. 

J J 

The arc (p,q) for the next iteration is obtained 
as follows. 


Let (u,v) be the leaving arc at its respective bound. 

If (u,v) = (i,0/ e W . set (p,q) = (r,s) 

If (u,v) = (i,j) e ( v «p q n W rs ^ u u rs )> (p,q) <ioes not char# 

If X pq = 0 or x pq - h pq , then set c pq - -b pq , ^ = -d pq , 
status = - status. Repeat this step. 

Step 8» The current basis structure is (T U i(p,q)], L,U). 

If c = 0 then the current solution X is optimum, stop. 


If c <0, then set c = - c , d = — d 

nn 9 r>n Th n * *nn 1 


pq 

status = -status. 


pq 


pq 7 pq 


pq’ 


If > 0, then set flag = 1 else set flag = -1, 

Go to Step 6. 

3.7 Numerical^ Exampl e s 

In this section a sm all-network flow problem is solved 
to illustrate the various steps of WC-lIF algorithm. The network 
is shown m Fig, 3«X. The numbers c , b d. are indicated 

XJ XJ XJ 

over each arc. -An amount of 5 units of flow is sent from source 
1 to sink 6. The aspiration levels and weights are specified 
as follows 

w ] _ = 0.15, w 2 = 0.10, r x = 0.35, r £ = 0.4, C ± = 165 and 
D 1 = 125. 
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The steps of the algorithm are summarized in Table 3.1. 
The ^ indicates the basic pivot arc and ^ indicates the basic 
arc leaving the basis. The graphs corresponding to the basis 
at various iterations are shown m Fig. 3.6. In these graphs s 
the basic pivot arc is drawn as dashed lines. The algorithm 
traced the path I, J, L, II as shown in Fig. 3.5 and obtained 
C * = 145 and D* = 155 as an optimum solution. 


3.8 Mathematical^ Formulation of_ Interval ^ Goal. Network^ Flow- 
Problem^ 


The mathematical formulation of interval goal network 
flow problem is as follows > 


Mm. 

s • t . 


Z = 


w 3 a 3 w 4«4 + r 1 P 1 + r 2 f 2 


2 x . - 2 x n _ 

( 0 ,i)el(i) J (i,a,eO(i) 1J 


(~ V > 

; v, 

Co, 


if i=l 

if i=n, V ieN 
otherwise 


0 

A 

X 

< b 


10 

1a 7 

z^ 

c . 

X. 

(i,3>e-H- 

U 

10 

Z 

d 

x- _ 

(i,0/£4 

10 

10 

Z 

c . 


(i,o) £ A 

10 

„ 

10 

Z 

d „ 

X 

(i,a)eA 

10 

10 


V ( l , 0 ) £ A 
+ ^1 ~ a l = ^1 
p 2 - a 2 = 

QCj ~ ^ ^ 

+ ^4 ~ a 4 = D 2 


(3.18) 

(3.19*) 

(3*2.0) 

( 3 . 21 ) 

(3.2^5 

(3.23) 

(3.24) 
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ig. 3-5 Paths traced by WGNF algorithm and IGNF 
alonthm for the numerical example 
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where, 

= lower limit on total cost. 

C 2 = upper limit on total cost. 

D-j_ = lower limit on total budget. 

D 2 = upper limit on total budget. 

We are interested in a solution between and C^. 
Therefore P^ and are minimized. Similarly if the solution 
is to be between and D then P 2 and are minimized. By 
assigning proper weights, we obtain the objective function (3.18). 

3.9 The_ Algor ithm. for^ jLnterval Goal_ Network Flow_ Problem* 

The algorithm for interval goal network flow problem is 
similar to that of weighted goal network flow problem except 
for a few modifications. 

The procedure to trace the path to reach an optimum 
solution is same as that of WGLF algorithm but the paths 
that these two algorithms take to obtain the optimum solutions 
are different. In interval goal problem, the decision maker 
specifies a range of aspiration levels for each of the goals 
and weighting factors for deviations from the ranges. The 
ranges specified form a rectangular region R in E whose four 
corner points are , (G 2 ,D^; , ( C 2 »D 2 ), (C^,D 2 ) as shown 

in Fig. 3.7. The algorithm obtains an optimum solution 
(X, c ty a^, p-p P 2 ) by minimizing the sum of weighted deviations. 




Fig, 3-7 Aspiration levels and devaiations for 
IGNF problem 



Fig 3 8 


Various paths traced 
to obtain optimum 


by IGNF alogorithm 
solutions 
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From above, it can be inferred that for different 
positions of R, the IGNF algorithm considers one of the 
corner co-ordinates of R as aspiration levels and obtains an 
optimum solution like the WGNF algorithm* Therefore, all the 
theorems proved for WGNF algorithm are also valid for IGNF 
algorithm. 

3*10 Nume ric al_ Example. 

In this section a small network flow problem is solved 
to illustrates the various steps of IGNF algorithm. The network 
is shown in Fig, 3.4. The aspiration level ranges and the 
weighting factors specified for this problem are as follows i 
C ± = 91, C 2 = 96, D 1 = 80, D 2 = 90 , w^ = 0.2, w 4 = 0.3, 
r-L = 0.15, r 2 = 0.35. 

The steps of the algorithm are summarized in Table 3*2. 
The indicates the basic pivot arc and ^ indicates the basic 
arc leaving the basis. The algorithm traces I,P,Q path as 
shown m Fig, 3.5 and obtains C =91 and D = 90 as an optimum 
solution. Q represents an optimum solution, in Fig. 3.5. 

3.11 Computational^ Results • 

The algorithms proposed m this chapter were coded, 
debugged and tested in Fortran -- 10 on DEC-1090 multi-prog ramm- 
ing, time sharing computer system. A number of randomly 
generated network problems were solved with different number 



Table 3.2 s Solution of the IGNF Problem 
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of nodes and arcs and computational times were noted. A 
network generator was used for computational study which 
generates well-structured networks. The program first 
generates a skeleton network for a specified width and lengtl 
and then adds arcs randomly until the network contains a 
specified number of arcs. 

Data structures based on augmented threaded index 
method[8] was implemented for storing the spanning tree. 

The tree was stored by means of thread indices, predecessor 
indices and the number of successors, each requiring an 
array of size n. Thread of a node i is a node that will be 
scanned in depth first search order of the sub-tree rooted at 
i. Predecessor of node i is the first node on the path from 
i to the root node of the tree. The number of successors of 
node i is the number of nodes in the subtree rooted at i. 

Link of a node i is the first arc on the path from i to the 
root node of the tree. If the arc is incident to i, it is 
stored as positive. If the arc is incident from i, it is 
stored as negative. 

When an arc leaves the basis, it results in the forma- 
tion of a hanging sub-tree. Thread indices help in scanning ! 
all the nodes of this sub-tree and thereby updating the dual 
variables. Predecessor indices are used for determining the 
minimum flow to be augmented, the leaving arc, and for updatii 
the flows in the cycle formed, when a non— basic arc is added 
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to the basis. The number of successors for each node is 
maintained. It helps in determining the common node for the 
two paths when traversed from two ends of the ente ring arc 
towards the root node. Once the common node is known, it is 
easy to determine the minimum flow to be augmented in the 
cycle, leaving arc and updating the flow as well as the basis 
structure . Link is an array of size n that stores the arcs 
belonging to the basic tree. The storage for preserving the 

j 

flows in arcs was reduce 1 ! from m to n. This was achieved by 
storing the flows of basic arcs requiring an array of size n. 

The remaining non-basic arcs exist either at their lower or 
upper bounds, respectively. If the non-basic arc is at its 
upper bound, its capacitor is made negative otherwise Its 
capacity remains positive. 

The mam emphasis with computational results was laid 
on (i) to check the number of iterations performed by the 
algorithm and (li) the computational time taken by the algorithms 
to get the optimum solutions. Problems sizes ranging from 10 
nodes and 40 arcs to 200 nodes and 2500 arcs which includes 
both the sparse as well as dense networks were considered for 
computational study. Arc capacities were randomly generated 
between 5 to 100 whereas c . ' s and d •' s were randomly gene- 

Xj X J 

rated between 1 to 50. Each problem was solved for three 
different aspiration levels and the computational times as 
well as the number of iterations are noted in Tables 3.3 and 
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3.4, About 20 problems of different sizes were solved for 
each of the algorithms. By comparing the columns (6), (8) 
and (l0) of Table 3.3, we can infer that the time taken for 
solving (C^, D^) e U is more than that for solving 

^ C l* D l^ £ ^ u r 4 or ( c i> D i) e s « $his is because, the 
optimum solution for (C^, D^) e U lies on B 2 U and 

thereby traversing a longer path. No comparison can be made 
between (6) and (8) of Table 3.3 as the time taken to solve 
(C^, Hj_) e U R^ and (C^, D^) e S depends on the values of 
and D^. One can infer similar results by comparing ooluians 
(6), (8) and (10) of Table 3.4 of IGNF problem. The reason 
being that WGNF algorithm and IGNF algorithm behave in a simi- 
lar manner. 

It is evident from the tables that both the algorithms 
can solve quite large problems in reasonable amount of time. 

For instance, problem of size 400 nodes and 2500 arcs was 
solved in about 2?r minutes by both the algorithms. Infact the 
time taken by both the algorithms are •omparable because IGNF 
algorithm is similar to WGNF algorithm except for a few- 
modifications which will not alter computational times sig- 
nificantly. 

Both the algorithms suggested by us are able to solve 
practically large problems in a reasonable amount of computer 
time. From the literature review we can conclude that no work 
was done In the field of goal programming techniques applied to 



Table* 3.3 s Computational Times 
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Table 3 »4« Computational Times of XGNF Algorithm. 

(Execution Times in Seconds on DEC-1090 Systems) 
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network flow problems. This made us to explore this field, 
and suggest exact algorithms for weighted goal and interval 
goal network flow problems. One disadvantage with goal 
programming techniques is, the decision maker should specify 
proper weighting factors otherwise he may not obtain the 
preferred solution. 
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KI 


go ra 270 

1 P ( : B h R . G E . 0 . OR * MUEM AX . LE . CBAR ) GO TO 230 

« 0 £«p*:BAR 

Tp |:3 ;| s u:-*|^.OR.<S|G«»|HARK< 5 E, 0 )GO to 230 

< " > ‘ - • ' ' . - ■ ' 

Wt?fc»s AXSMuE ' ■ . -:;• : : 

^ ‘ : - . - : . " ;■ . ' . 

Z . - ' R : ; 

S' ■ ■ rv^f * l > > i» 


m : , 

1 1 , w 

1 1 ti 

lW ' : 

;> ? j ; ■ 

:■! f • 

uj i >, 
j 3 5 s 

Lfi I b* 

1,1 IV 

I 1 » V 
1 1 !? , 

■ ■■ ■ 

i 1 2 » 

l!j? 
BUS 


r, : ' JSl > akio-, 

R • ' 1 2 s ) 2 V- 

FI.3 <* ( 1 )so 
FI , ?>:■- I S i K ) = V 
C = y * L A R 0 L 

1=51 MK 

9 0 K = f r‘ R ?;AlD ( I j 

PI ( iO = Pl ( PKf.O ( K) ) +COST(LINK(K)) 

P 1 0 ( K ) = P X IH P R f: L) ( K ) ) + 8 U 0 (1, 1 N K C K ) ) 

FLU'- {K ) 

1FCK.O0.P PKV)GO TO 100 

i — f\ 

GO n 9o 
100 Z 0 N f 1 F u f! 

CODE: =5 

110 ft s I s 'i = - ! 

CMI N 1*1 = 3 

0 1 C V I = o 

*'« ******* *+*♦*****♦♦* ************* ******** ********** 4 4 f * t 4 4 +»m 

1 * WaHWKiflK cw ** 1 " - ■ ■■> - 1 

*** »**<r'*«»*tf*t*** ♦*♦***♦#*♦**#♦♦*♦* ^4* ^ »<t },**$** + 

120 SIGM*-1 

«UE^AX=LARGE 

IFCDl .CjT.DMIMIM)5IGN = J 

130 FLO?<PQ = 0 • 

ITER*iTERt.l 

IFOl.EO.O.AMD.C.GE.CJ )GO TO 920 

1 F.l s 1 - & Q.C . AMD. 0. NE.D1 . ANO» (.CODE.E0 • 3 .CJR. CODE. EO. 4) J GO TO 
0 = ABS(.0-D1 j 

X F ( 1 GDDE.LE. 2 ) • AND, ( Q.LE.S) ) GO TO 500 

ENTERED 

DO 230 i = j ,;v 

STARrsPOIMT(I) 

FINlSH = POIrJT( 1 + 1 ) -1 
I F C 5 f ART . GT .FINISH) GO TO 230 
DO 230 JsSTART, FINISH 
K=H£AD( J) 

I F C GAP ( J 3 * LT . 0 ) GO TO 140 
CBAR*C0ST(J)+PI(I3-PICK) 

08A|*BU0 C J 3 + PIDC I ) •PID(K) 

STAlOSsi 

340 C R = b-PICD-COSTCJ) 

•* BMSS« 2 r?* PI0a) ’ BU0tJ) 

111 !|(:/]ak»Gf!o!5Hi ,g1.03G0 TO 230 

170 |fCJMR,LE.O,OR.OBAR.NE,0)GO TO 230 

C$J«*CBAR 
DMlMsDBAR 
5TA TPQ*STATUS 











? * ' - V ; ■ . ■•••••. •> s * * 

YHtS PAPr OF. THE.. PROGRAM DETERMINES THE LEAVING ARC AND THE 
AMOUNT .?•* FbJU i'Lf tit; AUGMENTED. 

* *• 1 *. *M:*M*M *** 4 *■* ***$***********$** + $ 4*4 

27 u iCA? = 7uP(r;uTERJ 

F L 0 P =UHSlCAP( PO) ) 

I F C S I* tl’Pa.EO. i )FEd«POsO 
l F ( F L. A 3 . E i> . 0 ) r; G TO 28 0 

F L 0 x p li = T F U 0 . v 

280 IFCICAP.LT.fDGO TO 290 
S T A i’ 7 6 2 1 
IX=i TAIL 
t Y-J HEAD 

DEL s* A 1 = IC A P-FLOWPO 
GO TO 300 
290 iCAPs-iCAP 
: ■; [ 

D M I ■ =-DMIM 
STAT USs-1 
IXsIHEAD 
I V = 1 TAIL 
DELTAlsFLOWPO 
3 • ■ ■ DELi A— 0 E L I A 1 

Rf LO/tfsf LOUPO 

310 I F ( 1 * „ £ Q . HO G 9 TO 390 

IFi J'iBERaxj .Gf. NUMBER U Y)) GO TO 350 
J = Li VK(IX) 

IF( J .{.1,0)50 TO 320 
CHA^3EsCAP(J)-FL0WCIX) 

30 T 3 3 30 

320 CHAN3£ 2 FLl)<*(IX) 

330 I F ( CHANGE . GE, DELTA) GO TO 340 

delta=:hamge 

K = IX 
LEG 2 1 

34= u^cm 

350 f§:f|^|J$J G0 T0 360 

BHA«3€*FLOvti(IY) 

3D n 370 

360 CHA«3€ 2 CAP(-J)-FL0W(IY) 

370 IF(C HAM GE.GE. DELTA) GO TO 380 

OELTAsSHANGE 

LEG 2 *! 

380 l?«PRCD(IT) 

so n 3iD 

390 ffiP n 

m&%ao TO 410 
STAi;a$=STATPO 

Z S3 1 pi I *ST ATUS*OELTA 
O'* IH 5 * 1N*STATUS* DELTA 
ID1H***J 

IFtOUGT.D)IDIR«l 
IFCOELTAL.lt. fCAPXSO TO 400 
GALL AUSMSN ; 

'‘mm. ■ 1 fJ 'i-% -4 4 A, . ■ ,v 1 

(50 j , * i u ■ 

4 00 Z' u'.. . HDiuTE .. A'.,./. . 

r* f> | m ; I 4 l ■ ■■■■ . . ... , 

•11 U ?F< DELTA j'.Lr.MAX)GO TO 460' ■' 

0 EDI ' ? = • ft X 

A,V®:- a o ‘ r, TAtUStOELTA 

Af.Dl =9 ^ S FA . US*OELTA 


b i. A l-i A ; 

x ! 

l./ i / 
F 2 J A' 
1(21 
p/i a. 
i 2 i o v 

bii- 


2ilr 











^ *• r > n ", 

U A K i :i K, 
f ; I'-'’ i ? b l; 



) = Q t ; { ) , / 



I ) .l.fi.S. AND. CODE. GE.3)G0 T3 420 

* ( w r L 1 , 0 f i j i f ,«l 1 f< I M , W i , W 2 / W 3 , N4 » 3 i/EH ( ID I R # D s, L < A , R i* t ■ A 0 i 
. 1 3 IHj ID 440 

tl ♦ M U X 


GO 13 920 

450 D p; L T A = Q E L T A 1 “ M A X 

4 b v hi 0=: M J ••' * S T A I i . S * 0 F: L T A 

A D D 1 = D M I K*ST AT US*0£LTA 
2 = C i ADD 
D = D t A DO I 
DELi DELTA I 

I F l DE L T A 1 . L E • S ) 3 0 TO 480 

470 CALL U PD ATI (C,Ct ,0,01 , hlf’] pi , W1 , *2 , W3 , W4 ,0 VER, IDI8, DELTA , RFLAG) 
C=C-A0D 
DsD-ADDl 

IFOVER.fiQ.l .AMD. t DELTA. EQ, ( DELTA 1. -MAX 3 3)30 TO 440 
IFt jVER.EO-.l 3Gu TO 92 0 
2 = 0 ADD 
D=D i ADD) 

480 IF(LEG.EQ.O)CALL AUGMEN 
IF ( LEG .ME .0 3 CALL UPDATE 
490 I~C3DE.E0.43GO TO 250 
GD T3 1 ?u 

¥4t{tft*t4.Mt*'M4 | t*****t**********t******** *********** *****'***** 

* THUS PAL ) J!< f’HliGP.M M SELECTS THE NON - BASIC ARC R5 FOR * 

* si mult;. " « i. ■ ■ mm (Aviol j >, * j u cycles, * 

************** ***4*****t********* ****** *******************'****** 

500 CBPQsCMIN 
DBP3=0M|N 
TAlLPQ=ITAIL 
HEAOPQ=IHEAD 

teO.-l )FLOWPO=-CAP{P0)~DELTA 

m m *-' . 


293b 


l>7 D 5 !. 
6 . 

t ?y7. 
tCC., 
v</9C 

its if 





WrmsWM- 

^w 

MM 


p; 0(1 ’/ ) 


. i ■ v . k : ■ ■ , , ■ . ■ . ■'■! J . I U i H . E y . - 1) U S' UA(, = -i 

i {• ( 3 P ’> , G . . - . . 0 , i L> X R , E y » 1 ) HFliAGs^t 

’i. P F f f,i A (J • n E . u ) & 1 G r* — « S LG N 
Fi,kG = ! 

r f ii j >-i - f b o ‘X p y 

GH Xj 7 . 4 { , 

I* kf (f. } » f ♦ I'l'i t ‘ t H- *t vl + r i 

■X THIS PA :•.;]• OF f H G PROGRAM DETERMINES THE LEAVING ARC A '■'•ID THE 

* ■ i S' F>.J • ro BE 'V Ki TED V* TWO Cl CUES SI * J»< ’ . '• * c '. 

4 4 ♦»<*»♦•♦*♦♦■»**♦***♦♦****###*♦»♦**♦*♦♦* **<'**•' 

5 -j 0 C D E 1 = i 

D K ij I ~ t* i i .) '» P y 

IF (Si M'PO.EO. 1 ) DELTA S =C AP ( PQ) -FLOWPQ 
LEG s o 

D- GT A2=IABSf CAP(PS) J/PATIO 
I F (DEL T A. 2 „ LT. DELTA1 ) D £ L T A 1 s D E L T A 2 

6 0 0 E M f £ R - P j 

S I AI'0S = STATPQ 
R I N = i 
INC*i 

I 1' A T L=T A I LPQ 
I HEX 0 = HEADPy 

610 I F t f ATOS.EQ.-l )G0 TO 620 
TX=: FAIL 
i'{ = i 16 AD 
Gn : o 630 
6 2 0 I X - i I E A D 

IY =) FAIL 

630 IPO K , EQ » I Y ) GO TO 790 

fp( HBER(IX) .GT.NlJMflER(IY) )G0 TO 710 
I F ( 2 0 0 E 1 .5 Q • 3 ) G 0 T 0 o 6 w 
J = LI MKUX) 

IF(J.{jT,0)G0 to 640 
SOMUX)=SirM(IX.) *RIM 
GO =’3 650 

640 S11K(IX)=SUMUK)-RIN 
650 IF ( CODE! • EQ, 1) GO TO 700 
66 0 IF (SUM (IX) 1670-700,680 

670 DGLTA2 = -FLOW(IX)/SUM(IX) 

GO I‘ 3 690 

680 J = lA3S(liIMK(IX) ) 

DEtms(CAP( J)-FLOW(IX))/SUM(IX) 

690 IFCDELTA2.GE.DELTA1) GO TO 700 


LEGSINC 
1C f ss X X 

TQT&LbFLO W ( K l ) -kDELTAI *SUM ( K 1 ) 
IXaPRED(lX) 

GO 5 3 630 

IFCODEi .GT.2JG3 TO 740 
Jssllil MK(Ilf) 

IF(J,LT*0)GO TO 720 
StlM(iy)sSUM(IY)-RIN 
GO 1 3 7 30 

SOM(IY}=SUM(IY)tRIM 
IFtCODEI «EQ. 1 3G0 TO 780 
XF( SUM <110)750,780, 760 % 
DEOTA2a*FLOWUY)/SUM( IY) 

{JO ,( 3 770 

|s 1 •% sS (L£ MX ( I Y) ) 

DEL TA?a(C4P( J 3-FLOW C I Y ) ) /SUM( IY ) 
IF< 0 €uTa2 .GE.DELTA1 )GO TO 780 
'hi SDELTA2 

.j^V .:«■«# m -*f |tj r*. 




■ ii 

' /llii : 


^ *Xi> v ,'i 


i r * < iu fi-iy 

)F»r 1 « 1 ) ; i/I’nH^aflD ' 

I r a I ii = T ft I l» B S ■ ■ 

I , J F ~\ D - - i t h j h ,H 

S I nl J f» s s f a D* s . 

F ! . U -V R > S F L :\ * n 5 + S T A frs+oelta ■ ■ 

i * " t ii t. 

E '' R -RS 
” ii L B ftOGRF'l 
D6 > . 1 >'. •= D SLi'AJ 
I I’ A! ij - S' A 1 9 j 
1 1-: t % k j ~>i E A D ?•,'/ 

r . . a - j* -t ; t m> 

F 1 0 1 - ■ P v ! = F l J r i P ,> + s r A T P o * D E I , T A 

I, t — >■ ft t 

EN J’ER=PO 
CALL A UG :-•£! ! 

IF (FLOORS. HE. 0. AND. FLOtoRS.he.IAB5(CAP(RS)))GO l*:i 320 
I F l A i S ( F L D«RS>. G E . S . A N D . A B S ( F L 0 u P S - 1 A R S ( C A P l k S ) 3 ) J G 0 T C 
S fAm>=-SrATPQ 
Cr PDs-CBPO 
D B P v’ = •* D 6 P 0 

I F ( r L 0 W PU • LE* S „ DR . ARS ( FLOWPO-I ABS (C API PO) ) ) „ LF . 5) GO TO 

S T A T P P = * S T A T P Q 

CAPCPP)=-CAPCP0J 

CBPi s-wBPP 

D B P ? = “ 0 8 P Q 

GD I'D 51 u 

F L 0 a P () = p b 0 vi P S 

CBP?=-GBRS 

DBP2*-DBRS 

PQ = RS 

TAlLPQftTAILRS 

headpo=headrs 

S I A! PQ*-STATRS 

IFCFLO^PQ.LE.O.QR.ABStFLOWpg-lAdSCCAPCPQ) )).LE.S)SO TO 

cap(pq)=-:apcpqi 

S1’ATPQ=-STATPQ 

CBPJJs-CBPQ 

dbp5=-dbpo 

GP 13 51 0 

IF(IABS(LEGKEQ.2)G0 TO 880 
ITAJLsTAILRS 
IHEA DaHEAORS 
STA!*JS=STATRS 
DELI’IAsDELTAi ♦RATIO 
FLOffRSa FLOORS* STAIRS ♦DELTA 

InTERsRS 
CALL AUGMEN 
ITAl L»f AILPQ 
IHEA9«HEA.OPO 
$r*rtJ5*5TATPQ 
OELTAiaDELTAl 
1 13*1 SI 
ENTERaPO 

c^imscbpo 

OMIMaDBPO 

IF(SfATUS*EQ.l )G0 TO 850 

5mI||s*5mJn 

ICAP*IAS5(CAPCPQ)} 

, G'l” ;r : -I . : .ILL 

RfLpfafliCWPQ 
n LU UPDATE ■' 

* , •« ... ( • >s-IAbS(CAP( LEAVE) 3 


I j; 


P ! i : 


' ** % ) li A «'**' iic 




5 

c * 

; , 
k . ■ 1 

C T; 

£ I. v/ ■, 

..Is,. % 


; j a a .'Pi; 


■ =0t:brAi^RA.no 

. a 2 




/ i 


5 5 "7 

5 !> . ■ 


8 9 u 




•*SS; 


"55 7 
, 5‘j / 

■ *> 5 1 , 
5 ■ 

i "5 1 

5 u . > 
,562.«, 
i>5f> 3 
>5646 
>5 t» 5 U 


960 


9 i U 
920 

930 


uti 

,! V. ■ 

( ; 9s LOG /2 
- ! “IR= 0 S 
K =t i 

CM = C B RS 
D * \ =0b»\S 
l‘F( 5 I a T tlS. 

c i 

D’ • =-D , 

3 C 5 P = J AnSf TAP(RS) ) 
t< £ il J •'< S F hi) )j HS 

c o-i updat* 

Z ■■ • . ) *• S(CAP( LEAVE) ) 

. . LE . S ) c AP ( LEA v E ) a-c AP ( LEAVE ) 

C • : ' • A)l.PU)-PI(HEA0P0) 

UP > , =BU0(P6) tPlD ( TAT LPV ) -p ID ( HEADED J 
i F J S 1 A'J’PO , S() . I 3 GO 10 90 u 
C0P/=-CBP6 
Of Pu*~OBPO 
S AIPQs-STATPQ 
Chf*u = ~CBP( 

DBP > = -DBP! 

3D O 810 
CAL! ('! J' l M E ( T I * 1 E 2 ) 
ni 13 ’ 5 (. 22,93 0 )C,D, TIME, JTEP 
FOR A n ! X , 'COSTS* “ ' 

s rop 

END 


El 4 . 4 , ' BUDGETS', FI 4 . 4 , 


TIME: 


I !i , ' T PEH = 


, 18) 


J-V 


¥ ¥ ¥ '¥ $ 4 1- 

* f 1 1 3 

* 4r‘lGH 
f f 4 I* ¥* t< 4 


: 7 


41 -■ 

' ‘ i f> , 

A+* 

JV ? * . 


W*- F 

* '? i 

"i M 4 I , 

a* ;/i,/ 

‘, 4 ;- I - 

y 8 2 S / 

„ v y 6 gji» 


io 

7 o 


8 0 
4 0 

30 

SO 


p«~ , - 

f% r r\ f ■ c* t *. * *.*• \ ....... - . ; ' ••• •...••/ •' 

* M- * * *■ t * '■ * t i *********$*#**** ******* v * ** y * i * * J: ♦ * * < M- 

StlBfOjritfE UPDATES THE OBJECTIVE FUNCTION OF FJKEU * 

'■' % ' .r j * 

♦ M'*« M»tM» *♦********************♦:** 4********** ******, ** 

S4;i-OT*’r ^ UPDATl (C,C1, 0,D1, MINIM, W1 ,W2 , N3 , W4 , 0-/tIR , 10IK,0* 
i RFLAu) 

REAL VIM, UN'IM 

m’sssrt over 

I *•' t 3 2 L i‘ A . Q E . Cl , 0 E * 6 ) ) 00 TO BO 
iF l *FuA3) 3 0 * 5 , 1 ) 
t F ( I 3 t R ) i 0 . 8 0 , 3 3 
tPCZ.sr.CUGJ TO 40 
4 f .*< J*(Cl-CJ M2*CD-D1 ) 
j 7 E R ~ 1 

ir t 'a j. /» r,Uri)Go to »o 

1 c 0 ]' *4= - 1 1, w ' ■■ 

3 /SR = 0 
RE TORN 

'.T.; = -n*(oci )*vmo-on 

:;n T j / r> 

l^C.GT.CUOO to 50 
ir.i-i 3*(Cl-C)tR4*(Ol~D) 

M f] ?0 

•i t < - «' ? f ( C *2 1 j f f) 4 ♦ ( 0 1 «*D ) 

so ro i<) 

£ 0 0 


i.tT’fi, 



IlSlllti 

HM 


r( !+i )-i 




J 

« 

6 t 

i 

» 

» , 

H 

y 

V ;# 


*■ * * * * f- * ********»*********, ***** 4 ********** 

I' ;. IIS Pf< -ii'i OH-UNS 4.1 OPTIMUM SOLUTION FOR |M f F;R VAL 13 
* P» ' f = 1 * p '> J H I -CL I TER T A NETWORK FLjw PROBLEM. 

* * '• ? •" •> = ■> i - " ■• • ’ ' ■ i ; • . - ,- 

I - 1 • : « - i I O 7 M, SEE D | ROOT, SOURCE# SINK. PREV , COUNT, F R 3 

1 F T i ISH, EN TER , V , STATUS , Q . PP , 1 , K , P , 

2 r i; m f; i , t i m e 2 , t i m e , C m l n , t , 

3 RUE.; ( 4 C. 5 ) , fA REAS) ( 4 05 ) , PI f 405 ) , PQ1 NT 1 40 5 ) , ? 

SCVK4-J5) 3 0 0 0 ) , C A p f 3 L 3 0 ) , C 0 S T ( 3 n 0 J ) , 

5 P M , u 5 , T A I L p 0 , > ! K A D P 0 , T A f L F< S , < i E A 0 R S , C 3 D K . C 0 0 


1 i V t V -> > I n.r, HU V 3 V V V ^ Hr < JU U } , ^ Ifl • l 3<<1) J ) , ) 

5 90, R S,TMLP0,'0F: A D P 0 , T A I L R S , H E A 0 R S , C 3 D F, , 2 3 3 E 

R >■: A L -1 A X , AX1.MAX2# M U E , M UKM4X , ICAP , I Fl, 0 «' , 3 1 3 I -1 , j F L 3 
0* • -ISIOR I i .• 1 fu': R ( 1 u S ) , L, 1 M-: ( 4 0 5 ) , f , f S I r A ( > 5 J # F ( , 0 ' ( -I 3 3 ) 
C 0 •“ 3 3 I r U l. , I H EAR , S r WHS , Di-XTA # XCAP.LEL , ENTER , C M I H 


1 LINK , FLO*. PPEU, SUM , CAP ,JHlMBEfC,'^<£Af>, P* JPI D F ' ,K/ 

2 , L E A v E F R F L 0 w 

rx pe iu 

10 FO R i A T ( 1 X » '>>*$) 

AC C 1' P r * , * I !>3‘!l , LENGTH , M , SEED , V , C 1 , 0 1 , C 2 , 0 2 
CALL 5ETRAM (SEED) 

CAUL US4( WIDTH, LENGTH, N,M, SOURCE, SINK, P01M T .HEAD) 

DO 23 3*1,3 

Z>.?\ J ) = IRA !(5, IDO) 

COST! J)=[RANC1 ,50) 

20 BOO ( J)= IRAN (1,53) 

Rl- A D t 2 i , *)*1 ,W2 f W3,W4 

CALL RTIMEITIMEI) 

* Tills PART 3F THE PROGRAM FORMS THE BASIS TREE ALONG iUTH f Ht: 

J PRFLjHCK.. i 1 DICES, HO. OF SUCCESSORS h vi) CAl-CULAfES * 

* THE DUAL VARIABLES PI'S AND PID'S. * 

OVER=0 
F L A S * 0 
S=1 .OE-5 
RFLA3=0 
f F L 3 (v * 0 
DO 3D 1*1,3 
PREDUJsO 
PI. (1)*0 
P1D(I)=0 

30 continue 

LARGS* » 000000 

”*0 

SlMI-«#1.0E + i4 
ft II D C 3 3 =0 
ROOTsSOURCE 
M A 3 M f 1 

CAPCWIsLARGE 

COST (MAO »LARGE 

ftflDtMAJaO.O 

PftEOC SINK ) *SO!JRCE 

LENKCSINK)sMA 

PR£tf»SX?<K>TOP«i t COUNTS© 

KOUNf si 

SC A X ( TOP)*ROOT 
40 IF U'OP, |0, 01 GO TO 70 

IaSC AN (TOP ) 
f F ( I •Gx 0) GO TO 50 

■ NUMlERU ) «COUNT*JLISTC I ) 

» f > ft m+ | » fi n ^ | . ; 

50 lii $ rc i i a-coyMT 

COyNfteCOUNTf 1 

AOCpRSV}*! 

op 3 *-l . .L'v 

yr.GE,N)GO TO 40 


. | i; u jj J w n r # i 

SUM , CAP , NU MRE R , 


I 

: > 4* v j. -h 

a; j* .. u 

' 9 ** 1 - - 


i 5 ( : i 1 1 l > , 

3 ![>{.•..* 3 '*• 

E 1 , COD!-: 2 , 1.) 
0« 

) , 5 1 5 M t A ; ) 

r ‘ . , 


<£ Ai) , P I , »1 3 , IF , K 



r^Kj=v 


I=SiMK 

90 K=ruR£AD(n 

Pi ('< ) =P X ( PRFD ( K ) ) ♦ COST ( DINK ( K ) ) 

PlD< <)=PI.>(PREiHK))+BU D ( L l K K ( K ) ) 

FIX v ( K i ~0 

IFfX.EO.PRElOGO TO ICO 
I = K 

GO n Pi) 

100 CONTINUE 
Z n 0 £ * 5 

1 10 RSI^ S **1 
DMINIHsD 

************ ********* ******** ********** ****************** ** ********* 

* ’ ft. ftft. V THE PKDGRA-i SELECTS A I" Ai>~ ■ ’ " > 'X C ~ > * 

* ( FROM^A^SEJ **♦*?♦♦♦**♦*♦♦*♦*♦*#♦♦♦♦♦*♦**♦**♦* - , * .. ,. ■ v I 

HnEMAX=LARGE 
Iff D2.LiT.DMINlM)SIGN = r 
ISO FLOaPQsO 

IF(9«INIM.GE.i:'l . AND.O«INlM.LF.D2,fr.NO.C.GE.CJ .AND, CODE, > GO 
1 88 0 

IF( ( ASS (D -1)1 ) .L.E,S.0R.ABS(02-D) .LE.S) .AND.C0D8 . LE . ? ) GO TO 4 60 
,, . S , rV-‘ . ( CODE. FQ. 3 * OR « C 1 ) J 1 rj «3i, 

I P C ( ABS ( 0-01 ) , LE . S. OR . ABS ( 0-02 5 . LE. S) .AND .(C.GE.Cl) 

1 .AND. (CODE, LE.2 ) )G0 TO 880 

IFfD.GE.D1 . AN D . D . LE . D2 . AND . C . GE . C 1 . AMD . C . UE . C2 ) GO TO 880 

ENTERsO 

ITER SITER+-1 

DO 230 1=1 ,N 

STAR TsPOlM T ( I ) 

FINISH=P01NT(I+1 )-l 
IF(START,ST.FINISH)GO TO 230 
DO 230 J=START, FINISH 
K'sHEADCJ) 

irC:AP(J),LT,0)GO TO 140 
CSAR*COST( J)+Pl (I)-PI(K) 

DBARS80DC J)+PIO(I J-PIDCK) 

STATUS'® I 

ft r 1 ') * i* ^ • - 

140 C8AR*PI (K J*PI ( 1 ) -COS'TC'J ) 

DARR*RIO(R)-PIO(I)-BUO( J) 

STATUS="1 
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